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Abstract 

We study a non-relativistic fermionic retarded Green's function by making 
use of a fermion on the Lifshitz geometry with critical exponent z = 2. With 
a natural boundary condition, respecting the symmetries of the model, the 
resultant retarded Green's function exhibits a number of interesting features 
including a flat band. We also study the finite temperature and finite chemical 
potential cases where the geometry is replaced by Lifshitz black hole solutions. 



1 Introduction 



At critical points, physics is usually described by a scale invariant model. Typically 
the scale invariance arises in the relativistic conformal group where we have 

t — > At, %i — > Xxi. (1.1) 

Here t is time and spatial directions of the space time. 

We note, however, that in many physical systems the critical points are governed 
by dynamical scalings in which space and time scale differently. In fact spatially 
isotropic scale invariance is characterized by the dynamical exponent, z, as follows [l] 

t -)■ X z t, Xi -)■ Xxi. (1.2) 

The corresponding critical points are known as Lifshitz fixed points. 

In light of AdS/CFT correspondence |2] it is natural to seek for gravity duals of 
Lifshitz fixed points. Indeed gravity descriptions of Lifshitz fixed points have been 
considered in |3](see also |4] for an earlier work on a geometry with the Lifshitz 



scaling,) where a metric invariant under the scaling (j 1 . 2p was introduced. The 
corresponding metric i^j] 

= i 2 (-% + fr + £'l • (i-3) 



where L is the radius of curvature^] The action of the scale transformation (1.2) on 
the metric is given by 

t -> X z t, Xi ->■ Xx h r ->■ A -1 r. (1.4) 

As a physical application, the Lifshitz geometry has been used to provide a 
possible holographic description for strange metals |5|^] In this setup the Lifshitz 
background is probed by D-branes with non-zero gauge fields in their world volume. 
By choosing the dynamical critical exponent, z, the authors of [5] have been able to 
match the non-Fermi liquid scalings, such as linear resistivity, observed in strange 
metal regimes. Having found the non-Fermi liquid scalings, it is natural to study the 
fermionic properties of the system to explore, for example, a possibility of having 
a Fermi surface in the model. To do so, one needs to consider a fermion on the 
Lifshitz geometry to find the retarded Green's function of the corresponding dual 
fermionic operator via AdS/CFT correspondence. 

Indeed utilizing fermions on asymptotically AdS geometries, it was shown that 
the AdS/CFT correspondence can holographically describe Fermi surfaces (7 13 



1 As it has been mentioned in jjjj, although the metric is nonsingular, it is not geodesically 
complete and in particular an infalling object into r = oo feels a large tidal force. 

2 In what follows we set L = 1. 

3 See also [6j for drag force computations in Lifshitz geometries. 
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Actually to see a Fermi surface, one should look for a sharp behavior in the fermionic 
retarded Green's function at finite momentum and small frequencies (for a review 
see e.g. |14|). Moreover the spectrum of quasi-particle excitations near the Fermi 
surface is governed by an emergent CFT corresponding to the AdS2 near horizon 



geometry of the black hole 10 



The aim of this article is to study fermions on the Lifshitz geometry which in 
turns can be used to study the fermionic retarded Green's function of the corre- 
sponding non-relativistic dual theory^ Fermions on the Lifshitz geometry or on a 
geometry with emergent Lifshitz geometry in its near horizon limit have been stud- 
ied in 17 and 18 respectively. In these papers, with a Lorentz symmetric boundary 
term, the real time Green's function of the fermion has been obtained. It was shown 
that the resultant Green's function has real self energy. We will come back to this 
point in section two. 

In this paper in order to find the retarded Green's function we consider the 
Lorentz symmetry breaking boundary condition introduced in [19[^J Since the cor- 
responding boundary condition preserves rotational and scale invariances, but breaks 
the boost, it is more natural to impose such a boundary condition on the geometries 
with Lifshitz isometry. Of course it also breaks the parity which is preserved by the 
Lifshitz symmetry^} 

Note that in this paper we consider the fermion as a probe. It would be interest- 
ing to extend this work to the case where the back reaction of the fermions is taken 
into account. 

The paper is organized as follows. In the next section we study fermions on 
the Lifshitz geometry where we find a solution for the equation of motion with a 
proper boundary condition. Then using the solution we calculate the corresponding 
retarded Green's function where we see that the model exhibits a flat band. In 
section three we extend our study to the finite temperature case where we show 
that although the system has excited zero energy fermionic modes at low momenta, 
at high momenta still it has a flat band. In section four we consider charged fermions 
probing a charged Lifshitz black hole where we show that while with the standard 
boundary condition the system exhibits a Fermi surface, in the non-standard case 
it still has flat band. The last section is devoted to discussions. 



2 Zero temperature 



The aim of this section is to study fermions on the Lifshitz background which will 
be used to find the retarded Green's function for the corresponding fermionic dual 
operator in the dual non-relativistic field theory. Before going into computations, 



15 



1G 



4 Fcrmions on Schrodinger spacetime has also been studied in 
5 Throughout this paper we refer to this boundary condition as non-standrad boundary condi- 
tion, while we refer to that introduced in [20j[2l] as standard boundary condition. 

6 The Lorentz symmetry breaking boundary condition has also been imposed for fermions with 



dipole coupling in 22 and, on the charged dilatonic black hole in 23 . 
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it is worth to note that the Lifshitz geometry is not a solution of the pure Einstein 
gravity with or without cosmological constant. 

In general to get the Lifshitz geometry one needs to couple the Einstein gravity 
to other fields. In particular the Lifshitz geometry may be obtained from gravity 
coupled to massive gauge fields. In the minimal case where we have only one massive 
gauge field, the corresponding action is given as follows 



1 



lGnG. 



where F, 



d„A u 



d+l 



d d+1 x yfg [R + A 



1 



\m 2 A 2 



(2.1) 



4 4 

It is easy to see that, with a suitable choice of the 



parameters m and A, the model admits the Lifshitz solution (1.3) with a non-zero 



gauge field given by 24 



A, 



2 z- 1) 1 



z 

,2 _ 



(2.2) 



For this solution the parameters m and A are m = Az, A = z 2 + (d-2)z + (d-l) 2 . 
Alternatively Lifshitz metric may also be obtained as a solution of the pure 



gravity modified by curvature squared terms 25 . As the simplest case consider a 



d + l dimensional gravitational action as follows 

/ = ^| d d ^xV=g(R-A + 



(3R 2 



(2.3) 



Using the equations of motion derived from the above action one can show that the 



Lifshitz geometry (1.3) is a solution of the equations of motion for a suitable choice 



of the cosmological constant A and the coupling constant (3 that are given by 

2z 2 + (d- l)(2z + d) 1 



A 



4A' 



(2.4) 



Although we could have Lifshitz metric in arbitrary dimensions, in what fol- 
lows we will consider the four dimensional Lifshitz geometry which could provide a 
holographic description for a three dimensional non-relativistic field theory. 



2.1 Fermions on Lifshitz geometry 

Let us consider a four dimensional Dirac fermion on the Lifshitz background whose 
action is 

S bu i k = J d 4 x^-gi^> l - (r a D a -DaV^-m * (2.5) 

where Jp = {e^ a Y^[d a + \{u pa ) a YP% with Y^ u = §[r", r"]. In our notation the space- 
time indices are denoted by a, b ■ ■ ■ , though the tangent space indices are labeled by 

fl,U---. 

Since the Lifshitz metric may be obtained from a gravity coupled to a massive 
gauge field, in general, the solution may also support a non-zero gauge field. There- 
fore one could consider a fermion that is charged under the background gauge field. 
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Nevertheless in what follows we will consider a neutral fermion. We will come back 
to charged fermions, latter, in the discussion section. 

To write the equation of motion one should use the variational principle which 
typically comes with a proper boundary condition. It is important to note that the 
boundary term is not unique and indeed there are several ways to make the varia- 
tional principle well defined using different boundary terms [19]. For the moment, 
we assume that there is a suitable boundary condition such that the variational 
principle will be well defined. With this assumption the equation of motion is 



(/ ,)'T"[^ + i(a;,„)„r^l 



m H> = 0, 



(2.6) 



where the nonzero components of vierbeins and spin connections for the Lifshitz 



metric (1.3) are 



(e t ) a = r z 5 ta , (e;) a = r5 l \ (e r ) a = r5 r 



and 



-Sin- 



Using these expressions the equation of motion reduces to 



F r z d t - ( - + 1 ) r r + rVdi + rY r d r - m 



* = 0. 



(2.7) 



To proceed it is useful to work in the momentum space where we may set 



,,iuit+ik.x 



ip(r). In this notation the equation of motion reads 



ir{k.T)ip 



-iur z Y l + (| + l) T r - rT r d r + m 



{2.1 



It is also useful to act by {Jp + m) on the first order equation of motion to find a 
second order differential equation which typically is easer to solve. Doing so, and 



using equation (2.8), one arrives at 



r 2 d 2 r -(z+ 2)rd r + u 2 r 2z + 



+i(z - l)ur z T r r + mV - m 2 



+ 1 



tp = 0. 



2Z2 



r k 



(2.9) 



In general this equation may not have analytic solutions. We note, however, that 
for a particular case of m = and z = 2 the equation has, indeed, an analytic 
solution. This is the case we will consider in this paper. In this case defining 
ijj± = |(1 ±FT*)^, one gets 



r 2 d 2 - Ard r + wV - r 2 (k 2 =f iu) + 6 



0. 



(2.10) 



4 



To solve the above equation we make the following change of variable 



ip ± (r) = r 3/2 e" r 7±(iwr 2 ) 



(2.11) 



by which equation (2.10) reduces to a well known differential equation for /±(£) 



d 2 f±(Q | d/±(fl | 



+ 



where 



A, 



± -. 

Aiu 4' 



e 



/±(0 = o, 



1 

4' 



(2.12) 



(2.13) 



We recognize the above equation as the hypergeometric differential equation whose 
solution is 

/±(0 = ct^±e-tF(a ± , -2 M± + 1,0 + 4^e^F(f3 ± , 2 M± + 1,0, (2.14) 

where F(a,b,£) is the confluent hypergeometric function, cf 2 are two constant 
spinors and 



± 



Therefore altogether we find 



(2.15) 



(2.16) 



± 



[IUJ H c 



4 Co 



with 

It is important to note that so far we have solved the second order differential 
equation and thus the constant spinors cf 2 are not independent and, indeed, re- 
stricting the above solution to be a solution of the first order equation of motion 
(2.8) leads to certain relations among them. More precisely one finds 



'lid 



-.T r (k.T)ci 



'IUJ 



-T r (k.T)cf. 



(2.17) 



We note, also, that the solution has not been uniquely fixed yet. In fact in 
the context of AdS / CFT correspondence one usually imposes a boundary condition 
at IR. In the Euclidean case the proper boundary condition is to assume that the 
wave function is finite at IR. When we are dealing with the real-time AdS/CFT 
correspondence, the proper boundary condition is to impose an ingoing boundary 
condition on the wave function at the horizon 26 . In our case using the asymptotic 
behavior of the hypergeometric function, 



F(a,b,0 



T(b) 



T(b - a) 



m e u a - b 

T(a) ^ ' 



for large |£|, 



(2.18) 
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the wave function is ingoing at "Lifshitz horizon" at r = oo, if the parameters 
and c\ satisfy the following relation 

by which the ingoing wave function near the Lifshitz horizon behaves as follows^] 

$ ~ r 2 e laj( *-^ + ^ lnr) . (2.20) 



One may wonder what "Lifshitz horizon" means! Actually the situation could 
be compared with that of fermions on the pure AdS case studied in, e.g. [27] where 
the ingoing boundary condition has been imposed at the AdS horizon at r = oo. 
In order to understand it better one may think of this boundary condition as a 
limiting procedure starting from a black hole solution and then approaching the 
zero temperature limit, as we will do in the next section. 



Alternatively to obtain the relation (2.19) and then the corresponding retarded 
Green's function^] one may use the prescription explored in |27j where the authors 
presented a derivation of the real-time AdS/CFT prescription as an analytic con- 
tinuation of the corresponding problem in the Euclidean signature. Indeed in our 
case, we have checked that, using this prescription we will arrive at the same results 
as those in this and the next subsections^ 

2.2 Retarded Green's function 

In this subsection we compute the retarded Green's function of a fermionic operator 
in the dual non-relativistic three dimensional field theory by making use of the 
solution we obtained in the previous section. One should note that, in the context 
of the AdS/CFT correspondence, in order to find the corresponding retarded Green's 
function it is crucial to appropriately identify the source and response of the dual 
operator. 

On the other hand the identification of the source and response depends on the 
boundary conditions which one imposes to get a well defined variational principle. 
Thus it is important to study the possible boundary terms one may add to the action 
to make the variational principle well defined. Therefore in what follows we will first 
find a proper boundary action for our model. To do so it is useful to explicitly fix 
our notation. 



7 Since in what follows we are interested in the low energy limit of the retarded Green's function, 
the momentum will be space like, i.e. k 2 < lu 2 . 

3 The prescription for calculating retarded Green's function in the context of AdS/CFT corre- 

33l. 



2<s 



spondence has been first considered in 26 and further studied in the literature in, e.g. 

9 Euclidean Green's function for fermions on the Lifshitz geometry has recently been studied 
in (17 . 
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Since we have been working in a basis in which PT* is diagonal, we use the 
following representation for four dimensional gamma matrices 



pr 



-a 2 



a 2 I' 
In this notation one has 



la 





ia 1 



i ' 



r 1 



-a" 








-a 



3 > 



-ia 2 

i 



1(7 



(2.21) 



1 



(i + rr')* = dia g (o, 1,1,0) 



*_ 



(l-PT*)* = diag(l, 0,0,1) 



/*A 




/o\ 


*2 




*2 


*3 




*3 


w 




W 






f*A 


*2 







*3 







v* 4 y 




W 



(2.22) 



(2.23) 



Therefore the boundary terms coming from the variation of the bulk action, 
5S hulk = l -J d 3 xV^h(^T r S^ - 5*r r *), 



(2.24) 



reads 



6S> 



bulk 



d xy-h 



-5*1*1 + 5*£* 2 + 5*3*3 - 5*4*4 



(2.25) 



Since the Dirac equation is a first order differential equation we are not allowed 
to impose the boundary condition on all components of the spinors. Thus the aim 
is to add a proper boundary term such that half of the degrees of freedom do not 
appear on the boundary. So we will have to fix only half of the spinors. 



We note, however, that the boundary terms may not be unique 19 . Of course 
different boundary terms lead to different physics. In our case since the dual theory is 
a non-relativistic field theory, one may relax the condition to have Lorentz symmetric 
boundary terms. 



Following the suggestion of 19 it is natural to consider the following boundary 
tern£3 

Sbdy = \ f ^xv^*!^ 2 * (2.26) 



which in our notation reads 



S U y = ^ / d 3 xV-h(^ 3 + *+* 4 - * 3 *l - *1* 5 



(2.27) 



'Note that this boundary term is different from that considered in 17 . 
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This boundary term is invariant under rotation and scaling but breaks the boost 
symmetry. Of course in our model, being Lifshitz geometry, the boost symmetry 
has already been broken by the geometry at first place. 

Adding this boundary term to the bulk action and varying the total action we 
arrive at 



5Sh u ik + 5Si 



bdy 



-h 



-h 



where 



(Xi,X 2 ) 
(Ci,G) 



1 

i 

71 



(*2 + ^4,^2 - *4)- 



(2.28) 



Therefore we get a well defined variational principle by setting a Dirichlet boundary 
condition on xi an d Cif^j As a result the source and response are given by (xi, C2) 
and (x2,Ci)> respectively. The retarded Green's function is essentially a matrix 
which maps the source to the response. To compute the corresponding retarded 
Green's function it is illustrative to explicitly write the solution we have found in 
the previous section in components 



W 



f D n F ( 



9 _ *w „2 

re 2 



<i;_; |; icur 2 



. 1. 

+ ' 2' 



/-> ' F (a . : -: /^/~ 
D+ F (a 



2' 



iujr 



\D^F (a_ : -: /^'/ 



. 1. 

-> 2' 



) + D n rF {p. 

) + K rF (£ 

) +J D+rF (ft 
) + DZrF (/?_ 



■ 3- 

-1 2 i 



+ ' 2' 



■ 3- 

' 2' 



■ 3. 

' 2' 



luor 



iujr 



(2.29) 



zu;r 



On the other hand in our basis one has 



r(fc.r) 



-h\(J 



(2.30) 



Therefore the equation (2.17) reads 



,+ 
-24. 
,+ 

-2t 



'%UJ 







/c 2 c n - A;ic u 




( \ 

Ac 



it 

V / 



(2.31) 



n If we had considered SSbuik — bSbdy the boundary condition should have been imposed on \2 
and Ci- 
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where the last equality is the ingoing condition (2.19) with 



= -2^^^. (2.32) 

By making use of this relation and utilizing the asymptotic behaviors of the solution, 
the retarded Green's function can be read as follows 



G R {k) = --r^—A . (2.33) 




As it is evident from the above expression for the retarded Green's function Gn(u, k 2 ) 
G 22 (u, —k 2 ) and det(G^) = —1. The spectral function is also given by 

A{k) = -hn(Tr(G R )) = im^l*!). (2.34) 
To study different features of the retarded Green's function one may use the rota- 



tional symmetry to set k\ — C In this case the retarded Green's function becomes 



diagonal. In fact taking into account that k = k 2 one gets 



'A-k 2 







A+k 2 

G R (k 2 ) = - I I . (2.35) 

£t 



Moreover one finds 



#.)= 2 4 g»!j l = (2.36) 



7T 



r 2 (^ + i)-r(r/)r(r/ + i); ' ' 



The behavior of the spectral function as a function of k 2 and co is shown in figure [T} 
Since the spectral function is symmetric under k — > —k, it is sufficient to draw the 
figure for positive k. As we observe the spectral function is positive for sign(w) > 
and diverges at u — > 0. Actually by making use of the asymptotic behavior of the 
Gamma functions one can read the asymptotic behavior of the spectral function 
near lo = 0. Indeed for finite ki one finds A ~ — showing that it has a simple pole. 

To further explore the physical content of the model it is useful to study the 
behavior of the eigenvalues of the retarded Green's function as we approach u = for 
fixed and finite k 2 - Indeed using the asymptotic behavior of the Gamma functions, 
for k\ 3> U), one finds 

k 2 - A kl k 2 + A to j . 

k 2 + A co k 2 — A ki 



12 Although the model has rotational symmetry, the resultant retarded Green's function (2.33) 
seems asymmetric with respect to exchanging k\ and &2. We note that it is the artificial of our 
asymmetric representation of the Gamma matrices. We will back to this point latter. 
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Figure 1: Three dimensional and density plots of the spectral function as a function 
of oo and k 2 . It is positive for sign(w) > and has a pole at u = for fixed k 2 - 




Figure 2: The real part of the flat band eigenvalue, Ai which shows a delta function 
behavior at u = indicating that the imaginary part of Ai has a pole at u — 0. 



This shows that for finite values of k 2 one of the eigenvalues, Aj, has a pole at u — 0. 
More generally one can see that the eigenvalue Ai has a pole at u> = for all values 
of spatial momenta. This can be seen, for example, from the behavior of the real 
part of the eigenvalue Ai where there is a delta function at u = as shown in figure 
[2] As a result one may conclude that there are localized non-propagating excitations 
in the model showing that the theory exhibits a flat band. 

Note that, as we have already mentioned in the introduction, the behavior of the 
retarded Green's function is different from that considered in 18 where by making 
use of a semi-holographic method it was shown that the self energy of the fermions 
appeared in the Green's function is real. We note, however, due to the fact we are 
using the Lorentz symmetry breaking boundary term the resultant Green's function 
has an imaginary part. Actually the situation is similar to the pure AdS case with 



Lorentz symmetry breaking boundary term (4.25) studied in 19 . Although in this 



case the bulk AdS geometry respects the Lorentz symmetry, the boundary term 
breaks this symmetry leading to a non-relativistic boundary theory. On the other 
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hand since the boundary term (4.25), up to parity, is invariant under the Lifshitz 
symmetry, the non-relativistic theory one gets from AdS bulk geometry has the same 
symmetry as if we had started with Lifshitz geometry in the bulk. Therefore one 
may conclude that the appearance of flat band is, indeed, the consequence of the 
non-relativistic feature of the dual theory. 



3 Finite temperature 

In this section we would like to redo the computations of the previous section 
for a non-relativistic theory at finite temperature. Following the general idea of 
gauge/gravity duality placing the dual theory at finite temperature corresponds to 
having a black hole in the bulk gravity. Therefore in our case we should look for a 
black hole solution in the asymptotic Lifshitz geometry. Actually black hole solu- 
tions in the asymptotic Lifshitz geometry have been studied in 34 37 . In particular 



the authors of 37 have analytically constructed a black hole which asymptotes to 
a vacuum Lifshitz solution with z = 2. The solution may be supported by different 
actions with different field contents, though the metric has the same form as follows 



ds 2 = -(1 - 



dt 2 



+ 



dr 2 



+ 



dx 2 



H 



(3.1) 



where r# is the radius of horizon. The Hawking temperature is 37 



T 



1 



2nr 2 H 



(3.2) 



3.1 Fermions on Lifshitz black hole 

Following our study in the previous section we will consider a neutral massless 



fermion on the Lifshitz black hole given by equation (3.1). For this geometry the 



non-zero components of vierbeins and spin connections are 



(e t ) a 



1 - 



(e,r = r,T, (e r ) a = rJl-^6 r 



H 



and 



(^tr)a — —{u T t)t 



1 



5: 



ta • 




'H/ V H 

Therefore the equation of motion for a massless fermion in this background, setting 



iuit+ik-x 2 



' x r ■ (f)(r), 



(3.3) 
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reduces to 



Or 



H 



- iurlrr 1 ) + iv r k ■ r 



o. 



(3.4) 



To solve this equation it is useful to act by Jp on the first order equation to get a 



second order differential equation. In fact, defining a new variable x 



dx 2 



dx 



1 -x 2 



one finds 



(3.5) 



where 



|(i ±r r r*) 



and 



+ i r H\ k 2 + uj 2 r\ 



Hi 



±- 



luor 



H ■ 



(3.6) 



The resultant differential equation has a well known form whose solutions are the 
associated Legendre functions P and Q. Therefore the most general solution of the 
above equation is 



±{r) = cf P(v,/j,±,x) + cf Q(v,n±,x) 



(3.7) 



where cf 2 are constant spinors. 

Of course so far we have solved the second order differential equation which its 
solution is not necessarily a solution of the equation of motion which is a first order 
differential equation. In other words the constant spinors cf 2 are not independent. 
In fact in order to find a solution of the equation of motion one needs to plug the 



solution (3.7) into the equation of motion which in general leads to certain relations 
between the constant spinors c 12 - Indeed using the recursion relations between the 



associated Legendre functions 38 , 39 , 



x~ 



dP(u, (JL±,X) 
dx 



+H±xP(v, fi±, x) 



-y/l — x 2 P(v,fi± + l,x) — jj±xP(iy } fj,±,x), (3A 



one finds 



c 12 = ir H T r (k ■ T)cf 2 - 



(3.9) 



In order to impose the ingoing boundary condition on the wave function at the 
horizon we note that at near horizon the oscillating part of the solution has the 
following form 



(3.10) 
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which in our notation the ingoing and outgoing waves correspond to plus and minus 
signs, respectively. On the other hand using the asymptotic behaviors of the associ- 
ated Legendre functions near x — 1 one observes that the function Q has both the 
ingoing and the outgoing components, though the function P has only the ingoing 
part. Therefore in order to have a physical solution one needs to set = 0. As a 
result the solution of the equation of motion of the massless fermions on the Lifshitz 
black hole satisfying the ingoing boundary condition is 

with c being constant two-component spinors satisfying 

c - =ir H T r (k-T)c + . (3.12) 



3.2 Retarded Green's function 



In this subsection, using the solution we just found, we will compute the retarded 
Green's function of a fermionic operator in the dual non-relativistic theory at finite 
temperature. As we mentioned in the previous section in order to compute the 
corresponding retarded Green's function one needs to properly identify the source 
and response of the dual operator which in turns depends on the boundary condition. 
In this section we will follow our notation in the previous section and will consider 



the same boundary action as that given by the equation (4.25). In this notation 



the source and the response of the dual operator are given by (xi,C2) and (x2,Ci)> 
respectively. 

In order to read the proper source and response one needs to find the asymptotic 
behavior of the solution as we approach the boundary. In fact by making use of the 
asymptotic behaviors of the associated Legendre functions (see for example [38]) one 

gets 



r 



l — V — fl± 



r(i + ^=p) 



r 2 ± ^iujt+ik-x 



A±r z c ± e 



2 ± Aut+ik-x 



(3.13) 



In other words one may write 



v|> ~ 



f A - c t\ 
A+c+ 

A+cf 



Auit+ik-x 



(3.14) 



fkicf + k 2 cf\ 



Note also that in our notation, equation (|3.12|) reads 

/ c t\ 




VI/ 



-th 






+ 



(3.15) 
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Altogether with this information the retarded Green's function of the dual fermionic 
operator in the finite temperature non-relativistic theory is 



G R (k) 



*j I k 2 A 2 _+2r H k 2 A- 



■k 2 A 2 _ 



2r H k 1 A-A + 
r 2 H k 2 A 2 _ 



-A 2 + 



V 



2r H kiA-A + 
r 2 H k 2 A 2 _-A 2 + 



A^+r 



%k 2 A 2 _ 



-2r H k 2 A_A + 



T k 2 A 2 



-A 2 + 



\ 



J 



(3.16) 



It follows from this expression that Gh(uj, k^) 
The spectral function is also given by 



(^22(^,-^2) and det(Gif) = -1. 



Aik) = Im(TrGfl) 

7T 



(r 2 k 2 A 2 + A 2 ,\ , x 

lm ( ^:_4 )- (3 - 17) 



To explore different features of the retarded Green's function it is useful to use the 
rotational symmetry to set k\ = 0. In this case the retarded Green's function reads 



G R (k 2 



/ r H k 2 A-+A + 
r H k 2 A_—A + 



\ 







th k 2 A- 



-r H k 2 A_ 



J 



Moreover for the spectral function one also finds 



(3.18) 



A{k 2 



-Im 



7T 



^r 2 a + x + )r 2 a + x-) + r 2 (i + x+)r 2 (i + x- 



r ^T 2 (\ + X+)T 2 (\+X-] 



r 2 (i + x+)r 2 (i + x- 



(3.19) 



± 



± rn\/k 2 + ui 2 r 



H) 



with X 

It is instructive to study the behavior of the spectral function in the small tem- 
perature limit. Physically small temperature means that we should look for the 



energies much higher than the temperature, 



i.e 



<C 1. Practically one may ex- 



pand the above expression for ior u 3> 1. Indeed by making use of the asymptotic 

,T 2 > 
^2, 



behaviors of the Gamma function, up to order of 0{^- 2 ), one arrives at 



A{k 2 



-Im 



7T 



X 



r r 2 (r / + l) + r(r / )r(r / + l) 

F»(i7 + ±)-rfo)rfo + i) 

zvrT (477 - l)r 2 (77 + |)r(77)r(l + 77)) 



Au r 4 (r? 



r 2 (r7)r 2 (77 + 1) 



(3.20) 



with rj — As we see at leading order it is exactly the same expression we have 
found for the zero temperature case (see the equation (2.36)). 

The spectral function as a function of uj and k 2 is depicted in figure [3j The plot 
is drawn for th = where the temperature is T = 1/4. To further explore the 
physical content of the model it is also illustrative to examine the behavior of the 
eigenvalues of the retarded Green's function as functions of u and k 2 . In fact the 
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Figure 3: Three dimensional and density plots of the spectral function as a function 
of u and k 2 at r H = h|=. It is positive for sign(w) > 0. Note there is no pole at 
oj = for finite hi. 



real and imaginary parts of the first eigenvalue of the retarded Green's function have 
been plotted in figure |4} As it is shown the imaginary part of the first eigenvalue 
has a pole at u — 0. This can also be seen from the delta function behavior of its 
real part. On the other hand it can be seen that the second eigenvalue has no pole 
at co = 0. 

In comparison with the zero temperature case we see that the spectral function 
has qualitatively the same shape, though there is a small deviation in the low mo- 
mentum modes. Nevertheless for high momenta it remains unchanged. Therefore 
the system has a flat band for high momenta. 

It is worth to note that at low momenta and for low energies there are several 
non-trivial peaks. In fact the presence of these peaks at low energies suggest that 
heating up the system has excited zero energy fermionic modes at low momenta. 



4 Non-zero chemical potential 

It is important to mention that when one studies fermionic features of a system 
in condensed matter physics, usually one looks for a possibility of having a Fermi 
surface. 

Actually in order to rigorously address this question one needs to consider a 



charged fermion propagating on a charged Lifshitz black hole^ where we could have 
a non-zero chemical potential. In fact as we have already mentioned in section two 
the Lifshitz geometry is not a solution of pure Einstein gravity. In order to find the 
Lifshitz solution one may couple gravity to a massive background gauge field. In 



13 Few days after submitting our paper, another paper 40 appeared on arXiv where charged 
fermions on the Lifshitz geometry were studied. Of course their background is different form what 
we are considering in this paper for charged fermions. 
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Figure 4: The imaginary (left) and real (right) parts of the flat band eigenvalue. As 
we see that there is a pole at u = in the imaginary part of the eigenvalue which is 
also evident from the delta function behavior of its real part. 



this case the background supports a non-zero gauge field. 

We note, however, that this gauge field diverges as we approach the boundary 
and thus cannot play the role of chemical potential. In fact in order to have a 
chemical potential, another gauge field is needed j4l]. Indeed the second gauge field 
has the proper near boundary behavior to define chemical potential. More precisely 



one may start with the following action 41 
1 



S 



d si 



-g 



R-2A--(d<j>) 2 



Ai4>/p(l)\2 



A 2 ^/p(2)\2 



(4.1) 



This model admits a charged Lifshitz black hole solution with critical exponent z 
for the particular values of A, Ai and A2 as follows 



A 



(* + !)(* + 2) 



Ai 



2L 2 ' -1 yfz~=V 
The corresponding black hole solution is 41~P 



A, 



Vz-1. 



(4.2) 



ds< 



2- t ?/2 1 dr 2 2 1 -2 
■ ' " + r ax , 



-r zz fdf + 



r 2 f 



f 



1 + rl {z+1) 

r z+2 



+ 



2(*+l) 

'_0 

r 2(z+l) ' 



K 



2(z-l) 



A 



(i) 



Azr, 



2(2+1) 




,2+2 



A 



(2) 



(2) 



(4.3) 
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Note that in comparison with the solution in 41 we have shifted the gauge field by constants 



to make sure that Af remains finite. Moreover by a proper rescaling we have set L = 1 
and also with a proper choice of the parameters the radius of horizon has also been set to one. 
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where vq and k are the only remaining free parameters which determine mass and 
charges of the solution, and 



i 



o 



It is, indeed, an asymptotically Lifshitz charged black hole whose Hawking temper- 
ature is 

2 + 2 A Z 2(z+l) 



T = 1 - -^ r ^ L> . (4.5) 

At low energy, using the general idea of AdS/CFT correspondence, the physics 
is governed by near horizon modes. In this case, for example, one should study 
fermions on the near horizon background. At zero temperature where Tq ^ = ^ L 
setting 

r-l = - — , t = -r, (4.6) 

(2 + 3z + z 2 )£' e ' 1 ; 

the near horizon background can be obtained in the limit of e — > where one finds 

ds*= s d : 2+d L, + d X i * 2 



(2 + 3z + ^ 2 )e 2 " 4(z + 2)' 

4(1) = A (2) = ^ (2) f4 7) 

T (z + lK' T (z 2 + 3z + 2)f ^ n 

As we observe the metric is in the form of AdS2 x R 2 . Therefore the low energy 
physics is described by an emergent IR CFT. Actually the situation is very similar 
to the relativistic case 10 . As a result we would expect that the model exhibits a 
Farmi surface whose physics is governed by an IR fixed point. 



4.1 Charged Fermions 

Let us first consider a four dimensional charged Dirac fermion on the Lifshitz back- 
ground whose action is 



S ] 



bulk 



1 



T a D a - D a T a )-m 



(4- 



Here Ip= (e^[d a + \{u pa ) a TP° - iqA^}, with = ^,T V }. 

As we discussed in the previous sections one needs to impose a proper boundary 
condition to get a well defined variational principle. With a suitable boundary 
condition the equation of motion is 



d a + 



m # = 0. 



(4.9) 
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The above equation of motion by the choice of ^ = (—h) 
to 



— 1/4 „— iut+ik.x 



~fV 2 r r d r - 



r zfl/2 



r + -r • fc 

r 



m 



ip(r) reduces 



ip(r) = 0. (4.10) 



As we already mentioned the low energy is governed by an emergent IR fixed 
point. To examine the low energy limit of the fermions one should consider the limit 
of u <C ■ At zero temperature, using the scaling (4.6) and in the limit of e — > 
keeping u/e fixed the above equation reads 



-m - i*T«(a; + f ) + *T • * - - , 

4 v 2 + 3z + z 2 



m 



m = o, 



where 



1 + z)^ 



e 



finite. 



(4.11) 



(4.12) 



2 + 3z + z 2 ' 

We recognize the above equation as a charged fermion probing the AdS2 x R 2 
background (4.7), as expected. As a result one may go through the construction 
of 10 to express the retarded Green's function of the fermion on the Lifshitz charged 
black hole in terms of the retarded Green's function of AdS2 model in small u limit. 

Here instead of doing so, one utilizes the numerical method to solve the equation 
of motion numerically. To proceed, it is useful to consider the following representa- 
tion for four dimensional gamma matrices 



— a'~ 








-a 



1(7 








lO 



1 ' 



r 1 



-a 








a 



-la- 



ta 



(4.13) 



Due to rotational symmetry in the spatial directions we may set &2 
using the notation 



0. Then 



$ 2 



(4.14) 



the equation of motion (|4.10|) reduces to the following decoupled equations 
1 



1 



r z jl/2 

for a = 1, 2. It is easy to see that 



® a ~ a a r r 



+ b Q r~ 



ia 2 + ma 3 



;-ir^y 



for r oo. 



$ Q = 0, (4.15) 



(4.16) 



To find the retarded Green's function, following 10 , it is useful to defined Ci = 
ipi/ip2 and (2 =i>3/i>4 where ^'s are defined via $1 = (^1, ip 2 ), $2 = (^3 5 V^)- These 
parameters satisfy the following equations 



r z jl/2 



ki 
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"f 1/2 d r (2 + 2m& 



n 



ki 



where 



r zfl/2 



Q = w + g/i (2) (l 



+ 



(4.17) 



(4.18) 



Using these equations the retarded Green's function is essentially given in terms of 
functions Gi(k,u) and G 2 (k,u) where 15 



G a (k,u) = lim r 2m C, 



for a = 1, 2. 



(4.19) 



with the ingoing boundary condition at the horizon which in our notation it is 1 10] 

Ca (horizon (4.20) 

The precise expression of the retarded Green's function in terms of G a depends 
on the boundary condition one imposes to get a well defined variational principle. 
For example, if we impose the standard boundary condition, in our notation, the 
corresponding retarded Green's function is 



G(k,u) = - 
Therefore the spectral function reads 



Gi(Jfe,w) 
G 2 (k,u] 



A(k,u) = -Im(Gx(k,u) + G 2 (k, 



UJ 



(4.21) 



(4.22) 



On the other hand for the boundary condition obtained by adding the boundary 
term (4.25) the corresponding retarded Green's function as a function of G a is given 
by (see also |22|) 



G{k,u) 



/ 2G 1 {k,w)G2{k,u) G 1 {k,u)-G 2 {k,uj) \ 

1 G 1 (fc,u>)+G 2 (fc,u>) G 1 (k,ui)+G 2 (k,u)) ' 

Gi(k,u))-G 2 (k,ui) -2 

\Gi(fc,w)+G2(fc,w) Gi(fc,w)+G a (fe,w) / 



(4.23) 



Thus the corresponding spectral function reads 

r G 1 (k,u)G 2 (k,u) 



A(k,uj) = 21m 



Gi(k,u) + G 2 (k,u) 



(4.24) 



Note that in the notation we are using in this section (see (4.13)) the boundary 



term (|4.25|) reads 



Sbdy — ; 
15 Here we set k\ = k. 



d z x^h(i)\i) A + ipli/; 3 + tpy 2 + ^i)(4.25) 
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Figure 5: The behavior of spectral functions with the standard (left) and non- 
standard (right) boundary conditions for z = 2 and T = as functions of u and 
k. For the standard boundary condition the system exhibits a Fermi surface at 
kf = 0.8902, while for the non-standard case it has a flat band. 



Therefore adding this boundary term to the action results to impose the boundary 
condition on a combination of the different components of the fermions as follows 

SS = i J d 3 x\^h(Sxlm - 8r}{xi - Xi^Vi + vt s X2), (4.26) 

where 

(Xl,X2) = ^|(^2 + ^4,^2-^4), 

071,772) = -^=(^1 + ^3,^1-^3), (4.27) 
4.2 Numerical results 

Having found expressions for the retarded Green's function and the spectral function 
for the cases of standard and non-standard boundary conditions, it is an easy task to 
find their behaviors as functions of k and u. Here to explore the physical content of 
the model we have plotted the spectral function of the model for both standard and 
non-standard boundary conditions. At zero temperature where f^ z = (z + 2)/z 
we set m = 0, qfx^ = \/3. For z = 2 the spectral function is shown in figure 
[5j While in the standard case the model has a Fermi surface at kf = 0.8902, in 
the non-standard case it exhibits a flat band. Note that since the retarded Green's 



function is an even function of k (see equations (4.15)) we only considered k > 0. 

It is worth to note that the model contains four free parameters which are mass 
m, critical exponent z, temperature T and chemical potential which always 
appears in the combination q/j^ where q is the charge of the fermion. Therefore it is 
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Figure 6: The behavior of spectral functions with the standard (left) and non- 
standard (right) boundary conditions for z = 3 and T = as functions of uj and k. 
For the standard boundary condition the system exhibits a Fermi surface , while for 
the non-standard case it has a flat band. 



natural to explore the physical content of the model when we vary these parameters. 
Actually changing the parameters we find the following behaviors. 

For fixed m, T and q^ 2 ' as we increase the critical exponent z for the standard 
boundary condition the sharp peak representing the Fermi surface becomes smaller 
and occurs at smaller kf and eventually for large enough z it distorts the Fermi 
surface completely. In other words the model does not have a Fermi surface. On 
the other hand for the non-standard boundary condition the model still exhibits a 
flat band though there is a depletion in the low momentum modes (see for example 
figure [6] for z = 3). 

For fixed z the dependence of the spectral function on the parameters m, T and 
g/z*- 2 ) for both standard and non-standard cases is the same as that for z — 1 where 



we have AdS black hole solutions (see for example 10 and 19 for standard and 
non-standard boundary conditions, respectively). 



5 Discussions 

In this paper, following the general idea of AdS/CFT correspondence, we have 
studied retarded Green's function of a fermionic operator in a three dimensional non- 
relativistic field theory by making use of a massless fermion on the asymptotically 
Lifshitz geometry. In this paper we have mainly considered the asymptotically 
Lifshitz geometry with critical exponent of z = 2. We have considered both neutral 
and charged fermions. 

Taking into account that the gravity on asymptotically Lifshitz backgrounds may 
provide holographic descriptions for strange metals [5] our studies might be useful 
to explore certain features of strange metals. 
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For neutral fermions at zero temperature where the bulk fermions propagate on 
the Lifshitz background the resultant retarded Green's function of the dual fermionic 
operator exhibits interesting behaviors. In particular we observe that the spectral 
function has a pole at u = for all values of spatial momenta. The appearance of 
the pole may also be seen from the behavior of the real and imaginary parts of the 
eigenvalues of the corresponding retarded Green's function. 

Having seen a pole at u = with all values of spatial momenta shows that 
there are localized non-propagating excitations which in turns indicates an infinite 
flat band. Actually the situation is similar to that in pure AdS geometry with the 



Lorentz breaking boundary condition 19 



We have also considered three dimensional non-relativistic theories at finite tem- 
perature. To study the three dimensional model at finite temperature we have 



utilized the asymptotically Lifshitz black hole obtained in 37 . We have shown that 
by heating up the dual theory, although the non-zero temperature can excite low 
momenta zero energy modes, at high momenta there is still an infinite flat band! 

An interesting feature we have seen in our model is that the spectral function is 
positive for sign(o;) > 0. In fact unitarity requires that the spectral function to be 
always positive. Since in our case the retarded Green's function changes its sign and 
indeed is negative for sign(u;) < 0, it is tempting to propose that retarded Green's 
function we have found for the non-relativistic model contains information for both 
particles and anti-particles! 

We have also considered charged fermions probing a charged Lifshitz black hole. 
While for the standard boundary condition the model exhibits a Fermi surface for 
non-standard boundary condition the model still has a fait band. We have also 
observed that as one increases the critical exponent z for the standard boundary 
condition the sharp peak representing the Fermi surface becomes smaller and oc- 
cures at smaller kf and eventually for large enough z it destroies the Fermi surface 
completely. In other words the model does not have a Fermi surface. On the other 
hand for the non-standard boundary condition the model still exhibits a flat band 
though there is a depletion in the low momentum modes. 

It is worth to note that in order to make the variational principle well defined 
one could also use another boundary action as follows 

Sbdy = I J d 3 XV^h (# + ^_-« + ) = % -j d 3 xV C ^(*J*l+*S*4-*I*2-*l*3). 

(5.1) 

With this boundary term the variation of the whole action leads to the following 
boundary terms 



-h 



Sxki + S(ki + x{8x2-(l5(2 , (5.2) 



where in this case the newly defined fields are given by 

(Xi,X 2 ) = -^(*i + * 2 ,*i-*2) 
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(Cl,C 2 ) = ^(*3 + *4,*3-*4). (5.3) 

If we follow the steps we went through in the previous sections one may compute 
the corresponding retarded Green's function in this case. Doing so, one finds that 
the resultant retarded Green's functions have the same form as those in the previous 
sections, except the fact that the roles of k\ and k 2 have been changed. Now in this 
case we could use the rotational symmetry to set k 2 = 0. Of course the physics 
remains unchanged after all. 

Note added: After submitting our paper we were informed by U. Gursoy that 
fermionic correlation functions on Lifshitz background has recently been studied 
in 42 . We note, however, that the authors of this paper have used different UV 



boundary condition than ours. 
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